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Abstract 



>- ■ The integr ability of two symplectic maps, that can be considered as discrete-time analogs of 

cn ■ 

Tjj- ' the Garnier and Neumann systems is established in the framework of the r-matrix approach, 



starting from their Lax representation. In contrast with the continuous case, the r-matrix 



Q ! for such discrete systems turns out to be of dynamical type; remarkably, the induced Poisson 

0^ ■ structure appears as a linear combination of compatible "more elementary" Poisson struc- 

tures. It is also shown that the Lax matrix naturally leads to define separation variables, 
JlLi' whose discrete and continuous dynamics is investigated. 

> 

X 



1. Introduction 

In a number of recent papers, it lias been sliown tiiat tlie well-establisiied teciiniques 
devised to construct integrable finite-dimensional hamiltonian systems out of integrable hi- 
erarchies of nonlinear evolution equations (often denoted "soliton hierarchies") [1] can also 
be applied in a discrete context, naturally leading to integrable maps [2-5]. As remarkable 
results, I would like to quote the construction of discrete-time versions of the Toda-lattice [6] 
and of the Calogero-Moser system [7] . 

In this context, the author has shown [5] that, starting from the so-called "Toda Hierarchy 
with Sources", one readily obtains integrable maps equipped with a Lax pair. In a further 
paper [8], the connection of these maps with the stationary Toda flows and with the flnite 
gap sector of the solution manifold of the Toda hierarchy has been rigorously established. 

In the present paper, we investigate the algebraic structure underlying two of such in- 
tegrable maps, namely the Discrete Garnier (DG) [3] and the Discrete Neumann (DN) [4] 
systems, in terms of the "classical r-matrix" formalism as generalized in a fundamental paper 
by Babelon and Viallet [9] . Surprisingly, we flnd that the associated r-matrices are of dynam- 
ical type, in contrast with those pertaining to the corresponding continuous systems [10,11], 
and moreover appear as a linear combination of "more elementary" r-matrices, so that the 
associated Poisson structure is itself a linear combination of compatible Poisson structures. 

We recall that r-matrices of dynamical type have recently been discovered for an ex- 
tremely important class of integrable finite- dimensional continuous time systems, namely the 
Calogero-Moser class [12]. 

A remarkable property of the maps under scrutiny is their separability in terms of "roots 
variables" naturally provided by the Lax representation, in full analogy with the continuous 
systems [11]. 

The DG and DN systems and their Lax representation are tersely reviewed in Section 2. 

In Section 3 we derive the associated r-matrices, while in Section 4 we define separation 
variables and investigate the corresponding dynamics. 

A few (possibly) interesting open problems are mentioned in Section 5. 

2. Two Integrable Maps and their Lax Representation 

The DG system is given by the following Lagrangean map: 

< ^„, ^„_i > ^„_i+ < ^„, ^„ > ^„+ < ^n, ^n+l > ^n+l = A^^ (2.1) 

In (2.1), n G Z, "^^ is an R^-vector of components \I/n , A is a diagonal matrix with distinct 
entries {Xi,- ■ ■ , Xn) and the symbol < ■, ■ > denotes the usual euclidean inner product in 



As explained in [5], the system (2.1) arises by taking A^ replicas of the Toda-lattice 
spectral problem [13]: 

a^.,^^\+h^^^^^+ar.^% = A,vl>W. (2.2) 

corresponding to N different values of the spectral parameter A, and by restricting the dy- 
namical variables (a^, bn) to the invariant manifold: 

a„ = 2< ^^,^^+1 >; br,=<-^n.-^n> (2.3) 

The corresponding Lagrange function reads: 

C{x, y) = < x,y >^ +- < XjX > — < X, Ax > 

where x = \l/n, y = "^n+i- 

Canonical variables q,p can be introduced in the standard way [2]: 

^ = ^n; P = ^\i^-^^-.y=^^) = ^^^^^^_^^ (2-4) 

thus allowing to recast (2.1) in the form of the following map: 

q' = {A- < q,q> q-p)/g (2.5a) 

p' = gq (2.56) 

with g'^ =< q^Aq> — <q,q>^ — <p,q>. The map (2.5) is symplectic for the standard 
symplectic form dq A dp. Exploiting the procedure outlined in [5], the map (2.5) can be 
represented in the following discrete Lax form: 

L = AL'A-^ (2.6a) 

where L, A are 2x2 matrices, given by: 



-(A/2+ <p,Rxq>) V< P^ q>{l+< P, R\P > / <P,q> 

-y/<p,q>{l+<q,Rxq>) X/2+<p,Rxq> 



(2.66) 



^ ^ I (A- < g, g >)/V< P^q> -V<p'^Q' >/V<P.q> ] /2 6c) 



where 

The meromorphic invariant function A(A) = det[L{X)] is the generating function of the 
conserved quantities: 

^ r 

A(A) = -AV4 + 5] -^ (2.7a) 



N 

A -A, 

ill 

Ah — At 



and an exphcit computation shows that they are mutuaUy in involution 

thus entaihng the complete integrability of the map (2.5). 

Similar considerations hold for the DN system [4], described by the map: 

^^^^^ +6,^n + ^^^^^ = A*n (2.8) 

where the discrete motion is constrained on the unit sphere 5"^: < '^m^n >= 1, and 
accordingly the Lagrange multiplier hn is determined to be: hn =< \1/ti, A\l/„ > —1. 
Eq.(2.8) is a Lagrangean map for the Lagrange function: 

£(x, y) = log < x,y > — < x, Kx > 
In terms of the hamiltonian variables g, p, defined as in (2.4) by : 



eq.(2.8) becomes the map: 

p' = g\\{A-b)g-p\\ ; ^' = ^[^^1^ (2.9a) 

||q|| = 1, <q,p>=- (2.96) 

which is symplectic for the Poisson brackets: 

{?, 9} = 0; {q,p} = I -q(g)q] {p,p}=pAq (2.9c) 
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It enjoys the Lax representation (2.6a), with: 

1/2+ <p,Rxq> - <p,RaP> /\\p\i X (2 10a) 

<q,RAq> /\\p\\ -1/2 - <p,Rxq>' ^' ' 

^^((A-P)/IIP|I -IM/IW) (2,10,) 

Again, the invariant function A(A) = det\L{X)\ is the generating function of the integrals of 
motion: 

^(^) = Ey^ (2.11a) 



A -A, 



with 



{P jqk-Pkqj ^^ 
Afc — Xj 






The conserved quantities Ij have been shown in [4] to be in involution for the Poisson brackets 
(2.9c). 

3. r-Matrix Formulation 

As it has been shown for the first time in [9] , whenever a hamiltonian system is associated 
with a Lax matrix whose eigenvalues are in involution for a given Poisson bracket, it enjoys 
an r-matrix representation of the form: 

{Li,L2} = [ri2,Li]-[r2i,L2] (3.1) 

We have used the standard notations: 

Li =L(g) 1 = ^Lie'(g)l 

i 
i 

ri2 = Yl nte' ® e'' r^x = ^ rk^e' ® e'' 

where {e^}f£i is a basis for the matrix Lie-algebra which L belongs to, and, in general, the 
coefficients rik will be functions on the phase space {dynamical r-matrix). 



For both the systems introduced in Sec. 2, the Lie algebra is obviously s/(2), and we will 
use the Cartan-Weil basis: 

3_/l 0\ + f0 1\ _ fOO 



- =lo -1 ■ " - 0^ "^"=1 ol (='•2'') 



rather than the Pauli basis: 

cr^ = a+ + a-; a^ = i{a- - a+); a^ (3.26) 

As a matter of fact, since L has an additional (rational) dependence upon the spectral 
parameter A, it has to be regarded as an element of the loop algebra Q = sl{2) (g) C(A, A~^), 
namely as a formal Laurent series in A with coefficients in s/(2): 

p 
L(A)= J2 ^-tA" (3.3) 

A:= — cxD 

Incidentally, we notice that for matrices (2.6b), (2.10a), the series (3.3) is actually uniformly 
convergent in any compact subset of the annulus maxj\Xj\ < |A| < oo. 
The trace form on Q, given by: 

(L) = res tr L(A) = tr L_i (3.4) 

allows to identify Q with its dual (the space of linear functions on Q), and to consider the 
r— matrix as an endomorphism R on Q, rather than as an element of Q ^Q. In such a "dual" 
picture, eq.(3.1) induces the following Poisson bracket between two functions on Q: 

{/,fi'}L = {L,[df,dg]R) (3.5a) 

with 

[X, Y]r= [X, R{Y)] + [R{Y), X] (X, Yeg) (3.56) 

and: 

R{X) = 5^r,fce«(e('=),X) (3.5c) 

As is well known [9], eq.(3.1) implies involutivity of the invariants of the matrix L. 
In the dual picture, it is convenient to look at the functions: 

fk = res^trX'^L^ = ^(L,A'=L) (3.6a) 
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From (3.5): 



{hJj} = midfk.dfjU) = (L,[i?(A'=L),A^L])+ {k^j)=0 (3.66) 



In our concrete cases, where L-matrices are given by (2.6b) and (2.10a), the functions fk are 
related to the invariants Ij (2. 7b), (2. lib) by the formula: 

We can summarize the previous result by the following Theorem: 

Theorem 1. 

The hamiltonian flows of the functions fk correspond to completely integrable continuous-time 
hamiltonian systems, endowed with the Lax representation: 

The proof of the above Theorem is straightforward. The complete integrability follows from 
(3.6b). As for the Lax representations, we have, by definition: 

8T 

a^:={L.A} (3.76) 

On the other hand, for any basis element cr^, and for any integer j, it holds: 



namely: 



and thus 



id{res X> L) 



tr a'^-— = tr a' {res A^[L, i?(A'=L)]) 

Otk 



BL 

resA^(— -[L,i?(A'=L)]) = 



whence the result, o 



It is worth to notice that there is an intimate relation between an integrable map and 
the hamiltonian fiows of its invariant functions; namely it holds the following Proposition: 

Proposition I 

An integrable map is a Backlund transformations for the hamiltonian flows of its invariants. 



The proof of the above assertion is trivial: let us denote by a; a point in the phase space 
M = (R^^jOj), and by x' = ^{x) a symplectic map possessing the A^ invariant functions J-'j, 
in involution with respect to the Poisson bracket induced by uj. 
Let: 



then: 



ux • ox 



but Tj is an invariant function for $, and thus: 



dx' 
and finally, since $ is symplectic: 



dt, 



{$(x),^,($(x))} 



i.e. $ maps solutions into solutions, o 



dx 

— = {x,j^j}\x>=cs>( 



x) 



We now proceed to the explicit evaluation of the r matrices. 

1) DG system. 

By direct calculation, we have: 

{L3(A),L3(;u)} = 



{L3(A),L±(^)} = ±^[L±(A) -L±(^)] 
A — /U 

{L±(A),L±(^)} = -^=L=[L±{X) - L±(^)] 

{L±(A),L^(^)} = ±-l-[L3(A) - L±(^)] + J_^ [L±{\) - L^{p)] (3.8) 

The above formulas already show the appearance of the "dynamical term" , ''" . They 
entail the following expression for the Poisson bracket {Li,L2}: 

{Li(A),L2(^)} = 
8 



^ 'U,L,{X)+LM] + ^ ,} =([a^0(a--a+),Ll(A)]-[((7--a+)0a^L2M]) 



(3.9) 
where 11 is the usual permutation operator: 



U = ^a'0a' 



Comparing with the general formula (3.1), we have the final expression for the r— matrix: 

ri2(A, //) = - Y^n + -— L=a3 ® (a" - a+) (3.10) 

To derive the endomorphism R{X) (3.5c) that takes part in the Poisson bracket (3.5a), first 
of all we recall [10] that the term ^^^H dualizes into the difference of the projectors V-^- and 
V- on positive and (strictly) negative A-powers: 



V+{X) = Y^XkX 



k 
/VfcA 

fc>0 



V-{X) = Y^XkX 

k<0 



k 



As for the dynamical term, we note that: 



resX ^L 



V<P^(1> 



-^<p,q> 

whence it follows: 

^<P.q> = resh-\L+ - L_) = +1(L, X-\a- - a+)) 
so that: 

R{X) =X.-X++ ^ TTT^3(a- - a+, X) (3.11) 

(L, A ^{a -a+)) 

It might be of some interest looking at the role of the dynamical part of the r— matrix in the 
Jacobi identity, that we write down both in the tensor picture and in the dual picture: 



[Li, [ri2,ri3] + [ri2,r23] + [r32,ri3]] + [i^i,{iv2,ri3}-{L3,ri2}] + cyclic perm. ^ (3.12a) 

{L,[X,B{Y,Z) + {{L,Y),R{Z)}-{{L,Z),R{Y)}]) + cyclic perm. = (3.126) 
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In (3.12a) all quantities are understood to belong to ^ ® ^ (S> ^ and, as usual, the subscript 
denotes the space which the corresponding tensor acts on nontrivially: 

Li = L (g) 1 (g) 1, ri2 = ^ Vjke^ (g) e'' (g) 1, etc. 

jk 

In (3.12b), we have shortly denoted by i3(-, ■) the usual Yang-Baxter term, namely: 

B{X,Y) = [R{X),R{Y)]~R{[X,Y]r) 
For the sake of simplicity, we focus our attention on eq.(3.12a), and write for a moment: 

(c) , (d) 

where r^2 stands for the constant part —^^^11 and r[2 stands for the dynamical part 

^^7^a3®(a--a+). 

We see the following: 

(i) As is well known r^'^' satisfies the Yang-Baxter equation: 

[ri2, ris] + [ri2, r23] + [^32, ^3] = 

(ii) The mixed terms, containing both r^'^' and r*-^-* vanish. 

(iii) The dynamical part r^'^^ yield both quadratic and cubic terms in -7= — =. The quadratic 
terms in the first commutator of (3.12a) cancel with the quadratic terms in the second com- 
mutator, while the cubic terms, appearing just in the second commutator, cancel among 
themselves due to cyclic permutations. 

We point out that property (ii) entails that the Jacobi identity equation (3.12a) splits into 
two equations, involving separately r^^' and r^'^\ which are both satisfied. In terms of Poisson 
brackets, this means that for the DG system the Poisson bracket (3.1) (or (3.5a)) is actually 
the sum of two compatible Poisson brackets, generated by A'^^ and r^'^^ respectively. 

(2) DN system 

The Poisson brackets between the elements of the Lax matrix now read: 

{L3(A),L3(^)} = 

{L3(A), L±(/i)} = ±-^[L±(A) - L±(^)] + ^L±(/i)(L±(A) - L^(A)); 
{L±(A), L±(^)} = ^[^±(A) - ^±(^)] + ^(L3(A)L±(^) - L3(^)L±(A)); 
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{L±(A),L^(^)} = 

= ±Y^[Ls{X) - LM] + ^[^±(^) - L^W] + T^[L^iX)LM -L± (^)L3(A)] (3.13) 
A — /U \\p\\ \\p\\ 

Comparing with (3.7a), we notice that hnear terms in L are the same, up to the substitution 
V^¥7q~> ^^ Iblh consequently, the r— matrix will be of the form: 

ri2 = --^ + -^(7^ (a- - a+) + r[f^ (3.14a) 

where r^2^ takes care of quadratic terms in L in (3.13). Skipping out the computational 
details, we report the result: 

rS?^(A,^) = ^(^+ + ^") ® W^Hf^)] (3.146) 

and correspondingly: 

rif\X.fi) = -^[aM{X)] ® (a+ + a-) (3.146) 

2\\p\\ 

It is also possible to write the dynamical term in an invariant form. In fact we have 

\\p\\ = —res tr La~ = —{L,a~) 
Consequently, the r— matrix (3.14) dualizes to: 

i?(X) = X_-X+-— i— a3(a--a+,X)-— i— (a+ + a-)([a3,L],X) (3.15) 

2(L, a ) 2[J^,o- ) 

The check of the Jacobi identity is now considerably more involved than in the DC case. 
However, the same remarkable phenomenon occurs : namely, the Jacobi identity equations 
decouples in three equations individually satisfied by r*^'^-', r^'^^'^ j-(d2)^ thus entailing that 
the Poisson structure characterizing DN is actually the sum of three compatible Poisson 
structures engendered by r'^-*, r*-^^-*, r^'^^'. It is perhaps worthwwhile to point out that the 
dynamical nature of r]^2^ ^ namely the presence of the ||p||~^ factor is indeed essential: the 
tensor [a^, L{X)] (S> {o'~^ + a~) alone is not an r— matrix. 



4. Separability 

In this Section we will show that the integrable maps denoted as DC and DN, as well 
as the integrable continuous flows (3.7b) associated with them enjoy the classical property of 
separability. Namely, we will show that there exists a set of canonical coordinates (/Uj, ttj)^^-^ 
such that: 
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where the function Wj depend just upon the variable iXj and the integrals of motion: 

Paraphrasing, with slight modifications, the derivation presented in [11] we will introduce the 
variables Hj as the zeroes of suitable polynomials taking part in the Lax matrix. 

To this aim we notice that the L_ element of the Lax matrices (2.6b), (2.10a) can be 
written as: 

where 

a = -^<p,q> (DG); a = \\p\\ (DN) (4.26) 

and: 

N 

g(A) = n(A-A,) (4.3a) 

P(A) being a monic polynomial of degree A^ for the DG system, of degree A^ — 1 for the DN 
system, whose real zeroes, which are functions of the canonical variables q, p, we denote by 

P(A) = n(A - h) (4-36) 

Following [11] we set: 

TT, = Lsifij) (4.4) 

i.e.: 

N 






-.-^ + E^^ j = l,---,iV-l (DN). 

Then, by direct calculation, starting from formulas (3.8), (3.9), and paraphrasing again [11] 
we can establish the following theorem: 

Theorem 2: The variables T^j^fij are canonically conjugated. Moreover, they are sepa- 
ration variables. 
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The separation equations are obviously given by: 



tt] = A(^,) (4.5) 

and consequently the functions Wj of formula (4.1) read: 

Wj= I ' dXVA{X) (4.6) 

We shall now look at the continuous flows of the invariants fk taking part in Theorem 1. 
From the Lax equation: 

87 

— = [L,R{X'L)] (4.7) 



we deduce: 



whence: 



Otk 



^^ |a=m. = -2Mi")(;U,)7r, (4.^ 



dtk 

It is easily seen that, both in the DG and in the DN case, only the constant part of the r- 
matrix contributes to M_ (Hr)- Indeed, a direct and simple calculation leads to the formula 
(here and in the following a superscript dot denotes differentiation: e.g. Q{Xk) = -^\x=xj'- 

Mi'=)(M=2«$:-^^ (4.9) 



taking into account that: 



we get from (4.^ 



dL- PiPr) dfir 



I A=/i 



a- 



dtk ' '^^ QM dt 



d^r ^^QiPr) Y- -^s ^(A 



E-Vt^ (4-10) 



dtk P{fir) ^l^r-XsQ{X 

Recalling the so-called Lagrange interpolation formula, which is a plane consequence of the 
residues theorem: 

we can write: 
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y — ^ — M = _yit^ (4.12) 

Eq.(4.12) can be immediately integrated, yielding to the Jacobi inversion problem (solvable 
in terms of Riemann functions) : 

y / d\ ^^ = -(t, - 4°)) y 4^— (4.13) 

^iMC) Q(A)v^(A) ''^Q{Xs) ' 

In particular, for the first nontrivial flows (/c = in the DG case and /c = 1 in the DN case) 
we get: 

y r''' dX ^ = -(t-t(°)) X (^^'--^ (^^)l (4 14) 

^7^.(0)) Q(A)v^(A) ^' ^""W^-iiDG)) ^^-^^^ 

The evaluation of the discrete-time evolution of the separation variables ^j is consider- 
ably more involved. The basic starting point is now the discrete Lax equation (2.6a), that 
entails both in the DG and in the DN case : 



to be of course complemented by: 



Ls{h) = -'^j (4-15a) 



Ls{p'j)=n'^ (4.146) 



Eq.(4.15a) allows one to express the quantities {p'j^q'j^} in terms of {a^j}? {"""j}; by inserting 
such expression into (4.15b), one gets a set of formulas relating {//^}, {tTj} to {/Uj}, {ttj} or, in 
other words, relating {n'j} to {Hj} through the map invariants. The explicit calculation are 
rather cumbersome but on the other hand straightforward, relying upon repeated applications 
of the residues theorem. Omitting the details, we just report the final result: 






^; - (M - A) = e - ^7^7' - 2 ^ ,, ,. ^X (4-16) 



where: 

e = 1 and j = 1 , ■ ■ ■ , A^ — 1 in the Neumann case 
e = and j = 1, ■ ■ ■ , A^ in the Garnier case 
M^Efc^fc ; A = EfcAfc 

By recalling that ttj = a/A((Uj) and getting rid of the square roots, one gets an algebraic 
equation of degree 2N{2{N — 1)) in the DG (DN) case, as the coefficients of powers 2N + 
2, 2N + 1{2N, 2N — 1) identically vanish. Of course, only half of the roots of that algebraic 
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equation will also solve (4.16), where a definite sign for the ttj has to be chosen. Once the Hj 
have been found, the frequencies Uj of the maps can be evaluated by hyperelliptic integrals. 
In fact, denoting by 9j the variable conjugated to /j, given by: 



7„ ^ J 7„ ^ 



^ k 






we have: 



1 /"Mfc 1 

z/,(/i,---,/^) = -J^ / TT^f^^/^ 

^V-^^fc (Aj-^)VA(^) 



and the maps linearize to: 



Ojin) = UjU + OjiO) 



Concluding Remarks 

We have shown that the r-matrix approach to integrable systems can be succesfuUy ap- 
plied to the discrete time case, starting from the Lax representation. It remains an open 
question whether one can construct different Lax pairs, possibly in terms of N x N matrices 
with a polynomial dependence upon the spectral parameter, leading to constant r-matrices. 
Another problem worth to be looked at is the role of the r-matrix as far as the discrete-time 
dynamics is concerned; in fact, since to an integrable map one can associate a (family of) 
interpolating hamiltonian fiow(s), engendered by a (family of) hamilton function(s) function- 
ally dependent on the invariants, one has to expect that the matrix A entering in the discrete 
Lax representation be expressible as well as a function of the matrices M^'^' that define the 
compatible continuous dynamics and are constructed through the r-matrix. Finally, a few 
words about quantisation : without going here into the stimulating problem of quantising 
symplectic maps, we just mention two points: i) once symmetrized with respect to p and q 
the classical invariants (2.7b) (2.11b) have a simple quantum version to commuting (formally) 
self-adjoint operators; ii) as illustrated in [11] in the continuous case, the separation equations 
can be quantised to a set of decoupled one-dimensional multiparametric ordinary differential 
equations of Schroedinger type. 
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